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Geodiversity, the natural range of geological, geomorphological and soil features, is thought to play a 
key role in the development of Australian ecosystems and evolution of their biota, due to the widespread 
occurrence of old soils with impoverished nutritional status and the comparatively restricted occurrence of 
fertile soils. While associations between soils and vegetation characteristics such as scleromorphy were first 
noted a century ago, modem theories propose evolutionary and ecological processes that shaped Australian 
flora and vegetation through interactions between soil nutrition, plant functional traits and flammability. 
Evidence in support of these generalisations comes mainly from site-specific empirical studies, surveys 
of plant traits and their associations, classification and mapping of land systems, analyses of regional 
environmental gradients and phylogenetic studies. The extent to which soils and the substrates from 
which they are derived place constraints on the climatic response of biota has important implications for 
understanding future responses of the biota to anthropogenic climate change. Yet, worldwide, there are few 
studies that examine the relative contributions of geological substrates and climate to variation in vegetation 
properties over extensive sub-continental regions. This study used a spatially explicit approach to examine 
the relationship between vegetation and geological substrates over New South Wales, a region of 80 million 
hectares in south-eastern Australia spanning a diverse range of geology, vegetation and climate. It aimed 
to assess the fidelity of major vegetation types to geological substrates and estimate the overall influence 
of geodiversity on vegetation composition relative to that attributable to climate. The spatial data were 
drawn from maps produced by geological and vegetation surveys. Geological maps were re-classified into 
16 broad units reflecting textural and mineral characteristics considered likely to be influential on plant 
growth. They represent a component of total geodiversity related to broad landscape-scale patterns in major 
bedrock and regolith types. Vegetation maps were re-classified in 16 broad formations reflecting structural, 
physiognomic and functional characteristics of vegetation, and into a larger number of 99 classes reflecting 
species composition. The two spatial data sets were analysed to determine the diversity of vegetation types 
within geological units and fidelity of each vegetation type to each geological unit. The relative influence of 
climatic co-variables was also examined using spatial surfaces spline-fitted to 30 years of weather station 
data. The results indicate a strong non-random relationship between vegetation and geodiversity, with most 
vegetation types restricted to a narrow range of geological substrates. Partial variance analyses indicated 
that the influence of geodiversity on vegetation composition was stronger than, and largely independent of 
the influence of climate. Consistent with current theories, sclerophyllous vegetation formations and classes 
showed a strong association with geological units characterised by low levels of mineral nutrients. It was 
concluded that landscape patterns of geodiversity are likely to place significant constraints on the response 
of native vegetation to future climate change. 
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RELATIONSHIPS BETWEEN GEODIVERSITY AND VEGETATION 


INTRODUCTION 

Gcodiversity is the natural range (diversity) of 
geological (rocks, minerals, fossils), geomorphological 
(land form, processes) and soil features. It includes 
their assemblages, relationships, properties, inter¬ 
pretations and systems (Gray 2004). Many of the 
earliest phytosociologieal studies in south-eastern 
Australia noted the association between plants and 
the geological substrates and land forms on which 
they grow (McLuckie & Petrie 1927, Pidgeon 1937, 
Fraser & Vickeiy 1939, Crocker 1944, Beadle 1948). 
The soils produced by weathering of these substrates 
in particular geomorphic settings vary greatly in levels 
and proportions of mineral nutrients, and in their 
textural and structural characteristics that govern their 
capacity to retain moisture and conduct subterranean 
oxygen. Different plant species vary in their ability to 
extract these three essential resources and to tolerate 
extreme levels of supply. Species that share similar 
ranges of tolerance to soil-related resources may 
therefore be expected to co-occur within communities 
on particular groups of geological substrates, at least 
in cases where there is vertical concurrence between 
soils and the underlying substrate. 

At biogeographic scales, geological substrates 
are thought to play a key role in the development of 
Australian ecosystems and evolution of their biota, 
due to the widespread occurrence of old soils with 
impoverished levels of nutrients and the comparatively 
restricted occurrence of fertile soils. Diels (1906) first 
noted a distinction between sclerophyll vegetation 
associated with sandy soils and ‘savanna’ vegetation 
associated with ‘favourable’ soil conditions. 
Andrews (1916) later proposed a connection between 
sclerophylly (a syndrome typified by small, thick 
leaves with thick cuticles and abundant sclerotic tissue) 
and soil nutrition, particularly nitrogen and calcium. 
By mid-century, empirical evidence was emerging 
that physiological and morphological traits related 
to uptake of phosphorus (Beadle 1953), acquisition 
of nitrogen (Hannon 1956) and accumulation of 
aluminium (Webb 1954) were closely associated with 
acidic, nutrient-deficient soils that occur in certain 
parts of the Australian continent. Modem theories 
propose evolutionary and ecological processes 
that shaped Australian flora and vegetation through 
interactions between soil nutrition, plant functional 
traits and flammability (Beadle 1966, Orions & 
Milewski 2007, Hopper 2009, Lambers et al. 2010). 
Evidence in support of these generalisations comes 
mainly from empirical site-specific or autecological 
studies (Beadle 1954, Webb 1954, Lamont 1982, 
Shane & Lambers 2005), surveys of plant traits and 


their associations (Gill 1975, Lambers ctal. 2010), the 
classification and mapping of land systems - areas of 
recurring patterns in topography, soils and vegetation 
(Christian 1952), analyses of local and regional 
environmental gradients (Mycrscough & Carolin 
1986, Keith & Sanders 1990) and phylogenetic 
studies (Johnson & Briggs 1975, Crisp et al. 2004). 

Sclerophyll elements of Australian vegetation 
appear to have originated on pockets of oligotrophic 
soils in the Cretaceous (Hill et al. 1999, Crisp et 
al. 2004), yet they did not diversify and rise to 
dominance until 25-10 million years ago. Their 
expansion and diversification was at the expense 
of mesic forest vegetation, and coincided with 
climatic cooling, drying and increased seasonality 
as separation of Australia and South America from 
Antarctica initiated circum-polar oceanic currents 
(Crisp et al. 2004). Expansion and diversification of 
the Australian arid flora occurred more recently, 5-2 
million years ago, as the continent moved still further 
north and experienced extreme wet-dry glacial cycles 
(Crisp et al. 2004). Fires also became prominent 
periodically and influential on vegetation during this 
period (Keeley & Rundel 2005). Yet strong edaphic 
patterns apparently persisted through this history 
and remain in the contemporary vegetation (Hopper 
2009) and, despite major extinctions and radiations 
associated with climatic upheavals, there is evidence 
ofbiome conservatism in a large majority of lineages 
(Crisp et al. 2009). 

While both soils and climate play prominent 
roles in theories of the evolutionary history of 
species within Australian vegetation (Crisp et al. 
2004, Orians & Milewski 2007, Hopper 2009), their 
historical inter-relationships are poorly understood. 
How much did soils constrain the historical responses 
of vegetation to climate change? This question seems 
crucial to understanding the future response of biota 
as anthropogenic climate change unfolds, yet the 
spatial dimensions of historical biomes are difficult 
to quantify and hitherto remain largely unexplored. 
Insights can possibly be gained by using spatially 
explicit data to study contemporary relationships 
between vegetation, soils and climate, yet this has 
not been examined at the level of assemblages over 
extensive bioregional scales. 

This study investigated the relative influence of 
geological substrates and climate on the biogeography 
of vegetation at a sub-continental spatial scale. Its 
aims were: i) to determine the fidelity between major 
vegetation types and geological substrates; and ii) to 
estimate the relative contribution of substrates and 
climatic variables to variation in species composition 
of the vegetation. The study was carried out across 
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New South Wales, a region of 80 million hectares in 
south-eastern Australia that currently encompasses 
a diverse range of vegetation, geology and climate. 
Historically, the region underwent profound climatic 
upheaval since the appearance of flowering plants 
in the Cretaceous. A spatially explicit approach was 
applied using map data synthesised from extensive 
vegetation and geological surveys carried out within 
the region (Keith 2004, Stewart et al. 2006). The 
two spatial data sets were first analysed to assess 
the fidelity between vegetation types and geological 
units. A second analysis was carried out to assess 
the relative influences of geological substrate and 
climatic factors on the species composition of 
vegetation. Both analyses were carried out on broad 
scale classifications reflecting on respective maps 
the broad structural, physiognomic and functional 
characteristics of vegetation and the broad textural 
and mineral characteristics of geological substrates. 

METHODS 

Vegetation map 

The vegetation map was assembled from 105 
source maps covering various subregions within New 
South Wales. These source maps employed a range 
of different vegetation classifications, varying spatial 
scales and overlapped with one another spatially to 
varying degrees. To simplify and standardise the 
source maps to a common format, the legend of 
each was re-classified into the vegetation formations 
and classes described by Keith (2004). In this 
classification 16 broad formations and sub-formations 
represent variation in structural, physiognomic and 
functional characteristics of vegetation (Table 1), and 
99 vegetation classes nested within them represent 
variation in vascular plant species composition 
(Keith 2004). Classes are not strictly nested within 
formations, as structure may vary considerably within 
compositionally defined units, however, classes were 
assigned to the vegetation formation representing the 
most commonly expressed structural, physiognomic 
and functional features of a mature stand (Keith 
2004). 

The 105 source vegetation maps were ranked 
according to their relative reliability based on 
assessments of classification skill, thematic and 
spatial resolution and currency using a protocol 
adapted from the one described by Keith & Simpson 
(2006). After standardising their spatial projections 
to the Australian Geodetic Datum 66 and Lamberts 
Conformal projection in a geographic information 
system, the source maps were merged sequentially 


so that features of the most reliable source map 
were displayed in preference to those of other maps 
wherever there were spatial overlaps (Keith 2004). 
A similar procedure was used to prepare a mask of 
extant native vegetation by reclassifying each legend 
category of each source map as native or non-native 
and then merging according to the currency of remote 
imagery from which each source map was derived. 
The mask was then applied to the composite map to 
derive an updated version (v3.0) of the vegetation 
map of NSW and the ACT prepared earlier by Keith 
(2004) (Fig. 1). 

Geological substrate map 

A map of geological substrates was derived 
from spatial data on the surface geology of Australia 
(Stewart et al. 2006) which, for New South Wales, 
was based primarily on mapping prepared by the 
Geological Survey of New South Wales (GSNSW, 
http://www.dpi.nsw.gov.au/minerals/geological/ 
geological-maps). This dataset shows the distribution 
of 785 geological units (mainly at Geological 
Formation level) and was generalised largely from the 
state digital geology dataset as at 2003, comprising 
the GSNSW 1:100 000 and 1:250 000 geological map 
series. Some areas in the north-east of the state (eg: 
Moree, Inverell, Tamworth, Manilla) and central-west 
(eg: Goulbum, Lake Cargellico) were re-compiled 
from more recent data sourced from GSNSW in 
2004-5 (Stewart et al. 2006). The basement geology 
of the Broken Hill region was compiled from 1996 
1:500 000 scale data from the national, NSW and 
South Australian geological surveys. Mapping for the 
Murray Basin region was compiled from 1991 1:1 
000 000 scale data from the Australian Geological 
Survey Organisation (AGSO). To compile a seamless 
state dataset, the original map data were edited 
along the edges of source datasets (edge-matching), 
which varied in age and spatial scale of compilation. 
Adjustments to some older geological data were 
made using geophysical data interpretation where 
particularly poor edge-matching or spatial accuracy 
(± 1 km) was identified in the source data (Stewart 
et al. 2006). 

The spatial data prepared by Stewart et al. (2006) 
generally did not distinguish contrasting depositional 
environments of unconsolidated Quartemary 
sediments along the coast. To resolve this deficiency, 
the coastal portion of Stewart’s map was overlain 
by more recent 1:25 000 scale mapping of coastal 
Quartemary geology, which discriminated a further 48 
units of sediments with alluvial, estuarine and coastal 
barrier depositional systems (Troedson & Hashimoto 
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Tabic 1. Vegetation formations of New South Wales. 

Formation 

Description 

Rainforests 

Forests of broad-leaved mesomorphic trees, with vines, ferns and palms. Includes 

Cunoniaccac, Sapindaceac, Monimiaccac, Lauraccae, Meliaecac, Myrtaccac, Apocynaccac, 
Rubiaceac, Aspleniaccac, Dryoptcridaceac. Coast and tablelands in mcsic sites on fertile 
soils. 

Wet sclerophyll 
forests (shrubby 
subformation) 

Tall forests of sclcromorphic trees (typically cucalypts) with dense understorics of 
mesomorphic shrubs, ferns and forbs. Includes Myrtaccac, Rubiaceac, Cunoniaccac, 
Dryoptcridaceac, Blcchnaccac, Astcraceae. Relatively fertile soils in high rainfall parts of 
coast and tablelands. 

Wet sclerophyll 
forests (grassy 
subformation) 

Tall forests of sclcromorphic trees (typically eucalypts), with grassy understorics and sparse 
strata of mesomorphic shrubs. Includes Myrtaceae, Poaccae, Euphorbiaceae, Fabaceae, 
Casuarinaceac and Asteraceac. Coast and tablelands in high rainfall regions on relatively 
fertile soils. 

Grassy woodlands 

Woodlands of scleromorphic trees (typically eucalypts), with understories of grasses and 
forbs and sparse shrubs. Includes Myrtaceae, Poaceae, Asteraceae, Epacridaceae and 

Pittosporaccae. Rolling terrain with fertile soils and moderate rainfall on the coast, tablelands 
and western slopes. 

Grasslands 

Closed tussock grasslands with a variable compliment of forbs. Includes Poaceae, 

Asteraceae, Fabaceae, Geraniaceae and Chenopodiaceae. Fertile soils of the maritime zone, 
tablelands and western floodplains. 

Dry sclerophyll 
forests 
(shrub/grass 
subformation) 

Dry sclerophyll 
forests (shrubby 
subformation) 

Forests of scleromorphic trees (typically eucalypts), with mixed semi-scleromorphic shrub 
and grass understories. Includes Myrtaceae, Poaceae, Asteraceae, Ericaceae, Dilleniaceae 
and Fabaceae. Moderately fertile soils in moderate rainfall areas of the coast, tablelands and 
western slopes. 

Low forests and woodlands of scleromorphic trees (typically eucalypts), with understories of 
scleromorphic shrubs and sparse groundcover. Includes Myrtaceae, Proteaceae, Ericaceae, 

Fabaceae and Cyperaceae. Regions receiving high to moderate rainfall on the coast, 
tablelands and western slopes. 

Heathlands 

Dense to open shrublands of small-leaved scleromorphic shrubs and sedges. Includes 

Proteaceae, Fabaceae, Myrtaceae, Casuarinaceae and Cyperaceae. High rainfall regions of 
the coast and tablelands on infertile soils, often in exposed topographic positions. 

Alpine complex 

Mosaics of herbfields. grasslands and shrublands. Includes Ericaceae, Asteraceae, 

Gentianaceae, Ranunculaceae, Poaceae and Cyperaceae. High, snow-prone parts of the 


southern ranges. 

Freshwater wetlands Wet shrublands or sedgelands, usually with a dense groundcover of graminoids. Includes 
Cyperaceae, Restionaceae, Juncaceae, Haloragaceae, Polygonaceae, Ranunculaceae and 


Forested wetlands 

Myrtaceae. Throughout NSW on peaty, gleyed or periodically inundated soils with impeded 
drainage. 

Forests of scleromorphic trees (eucalypts, paperbarks, casuarinas) with sparse shrub strata 
and continuous groundcover of hydrophilous graminoids and forbs. Includes Myrtaceae, 
Cyperaceae, Ranunculaceae, Blechnaceae, Poaceae. Floodprone plains and riparian zones 
principally along the coast and inand rivers. 

Saline wetlands 

Low forests, shrublands and herbfields of mangroves, succulent shrubs or marine herbs. 

Includes Verbenaceae, Chenopodiaceae, Juncaceae and Poaceae. Coastal estuaries and saline 
sites of the western plains. 

Semi-arid 

Open woodlands of scleromorphic trees (eucalypts, acacias), with open understories 


woodlands (grassy mostly of chenopod shrubs, usually with strong representation of perennial and ephemeral 
subformation) grasses and forbs, including many ephemeral species. Includes Myrtaceae, Fabaceae, 


Semi-arid 

woodlands 

(shrubby 

subformation) 

Chenopodiaceae, Asteraceae, Poaceae and Polygonaceae. Low-moderate rainfall regions of 
the near western plains on clay soils, including infrequently flood-prone sites. 

Open woodlands of scleromorphic trees (eucalypts, acacias, casuarinas), with open 
understories of xeromorphic shrubs, grasses and forbs, including many ephemeral species. 

Includes Myrtaceae, Cupressaceae, Fabaceae, Myoporaceae, Sapindaceae, Asteraceae, 

Poaceae and Acanthaceae. Low-moderate rainfall regions of the near western plains, 
including infrequently flood-prone sites. 

Arid shrublands 
(chenopod 
subformation) 

Arid shrublands 
(acacia 

subformation) 

Open shrublands of chenopod shrubs, with perennial tussock grasses and ephemeral herbs 
and grasses. Includes Chenopodiaceae, Asteraceae, Aizoaceae, Fabaceae and Poaceae. Low 
rainfall regions of the far western plains. 

Open shrublands of xeromorphic shrubs, hummock or tussock grasses and ephemeral herbs 
and grasses. Includes Fabaceae, Proteaeceae, Myoporaceae, Asteraceae, Casuarinaceae and 

Poaceae. Sandy or rocky landscapes in low rainfall regions of the far north-western plains. 
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Vegetation Formations 


Rainforests 

Wet sclerophyll forests (Grassy subformation) 
□ Wet sclerophyll forests (Shrubby subformation) 
Grassy woodlands 
Grasslands 


Alpine complex 
H Forested wetlands 
Hi Freshwater wetlands 
□ Saline wetlands 

]j Semi-arid woodlands (Grassy subformation) 


] Dry sclerophyll forests (Shrub/grass subformation) ( | Semi-arid woodlands (Shrubby subformation) 

Dry sclerophyll forests (Shrubby subformation) |H Arid shmblands (Acacia subformation) 
Heathlands IHI Arid shmblands (Chenopod subformation) 


Figure 1. Vegetation map of New South Wales (version 3.0) showing the distribution of 16 formations 
and subformations (see Table 1 and Keith 2004 for description) prior to intersection with extant mask. 


2008). In total, the combined spatial dataset mapped 
the distribution of 833 geological units across New 
South Wales and was re-projected to the Australian 
Geodetic Datum 66 and Lamberts Conformal 
projection in a geographic information system. As for 
the vegetation data, the map legend was simplified by 


assigning each lithological unit to one of 16 substrate 
types (Table 2). Substrate types were defined on the 
basis of attributes that were considered important 
influences on the supply of plant resources: mineral 
composition; weathering characteristics; texture and 
depth of derivative soils. The resulting map is shown 
in Fig. 2. 
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Table 2. Characteristics of geological substrate types in New South Wales 
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Low quartz Sedimentary & Mixture of clay minerals (e.g. montmorillinite, Clay & loam mostly Siltstones, mudstones, shales, Steep-flat terrain 

sedimentary metamorphic illitie, kaolinite, chlorite) mixed with silt >lm greywacke, phyllites, schists. 

particles of quartz, calcite, etc. producing soils Apsley-Macleay. 

with moderate-high levels of phosporus & 
exchangeable cations 


D.A. KEITH 



co 

CD 

6 

o 

£ 

o 


03 

fc 

CD 

03 

V 

Dh 

<D 44 
<d o 

55 S 


<D 

.p 

« s 

3 & 

<x> 

O in 
T3 


~ a 

o > 

a 5 

o 5 

4-* 03 

‘g £ 


CD 


p 

<D _ 

& p 43 
o p <d 

C/) XI rT) 


<3 

- a 

&. ,5 

£* A 
CO 

c c/5 

■B £ 

CL) 03 

*£ £ 
ta T 

CD 

o 

n ■a 

a P 
■§ 

2 go 
CO 

i 03 

(D 

C 


O 

o 


CO 

p 


.81.8 

a *-• cl 

£ ^ o 

43 <X> 43 
c.ti a 

C <4_( 

0.2 o 

^o 

co of ' 

2 © 

CD 43 

.£ £ 

S & 

O 

1 .s 

co O 

« o 

js§ 

4 =: o 
00 p 

2 a 


CD 
> 
£ 
43 

co 

_ 2 
cd 


£ .2 
■p 'co 
< CD 

£ cb 

• -H c3 


8 B 

'5h m 

^ 1 . 


xj 

P 

o3 

00 


a, 

(D 

CD 

X 

CD 


p 

CD 

'£ 

p 

a 

J-h 03 
CD ^ 
C Ph 
.53 co 

E ^ 

CO W 

O >, 
£ -o 

■ 2 Id 

£ P 
O g 

CD •£ 
(D 

< CD 

O *5 
>-l £ 


£ 

A 


c3 

£ g 
g; cn 
£* A 


c 

03 

CO 

2 

O 

43 

a. 

CO 

O 

43 

CL 


co co 
2 ^ 


J s 


p 

CD 


p 

p 

+-* 

CO 

O 

£ 


£ 

o 


> .2 
2 t3 

CD ° 

£ ia ~ 

X (D 0 

1 p 2 

- 03 <D 


O 


£ 

A 


a 

03 

co 

£ 

o 

43 

CL 

co 


CO co 

~ O 


CD 

> .s 

CD ' ' 


43 03 

CD ° 

£ - 

.§■§ 
£ o 

2 

£ c 

1 s 
> 43 

o >< 

h4 CD 


03 

•2 P 
cl 4 : 

£* v 

CO 




CD 

# n 

a 



> 

CD 

p 



03 

’S 




bb 

<D 






o3 

CO 

CO 

4-* 

CO 

£ 

CD 

O 

"P 

a 

p 

Q 

'co 

c3 

O 

.£ 

*£ 

t-i 

P 

o3 

co 


CO 

>, 

00 

^1 

o3 




CD 

£ 

P 

P 

c 

Id 


o3 

a 



U 

C/3 

00 

E 

0 

4-» 


43 

■ SP 
£ ^ 
.2 b 
■o £ 
c 

« JS 
to -33 

"2 S 

•§ 'o ’ 

a co 

o sp. 


CO 

« .8 
0 g 

D 

I'S 

J= ^ 
biD'P 
^ P 
03 


O 


O g- 
P3 .£ 
^ o 
. P 
■ X 

1 o 


03 

£ , 
C+H 

O 

co 

’d 

> 

CD 


43 t 

.SPX 

^ CO 

frig 

(D 

> ^ 


CO 


p 

O 

CD 

P 

bO 




co w 

•I I 
« « « 
2 

2 2 .s 


43 



43 ^ 

E -2 4 ^ 

4 —• QD 
' 1—1 P 


’> 

p 


■»-i > o3 

O *c 

CO 03 

CD C 

0 

co 

CD 

0 s 

co^ 

CD 

O E ’> 

CO rX ^3 
CD 9 43 

oi vE 

£ 

w 

Pi 

pi E.'p 


co 

S § 

• p o 

p .a 

s ^ 


co 

.£ 

’£ 

8 o. 
3 "g 
^ s 

O co 


03 


X3 

P 

03 

CO 


03 

CO 

c3 o3 

P 03 

P 03 

P 'X 

^ > 

S > 

£ P 

8 4 

8 8 

CO ^ 
CD P3 

Pi "3 

Pi ~3 

Pi O 


-p 

p 

03 


03 > 


O 

s 

. & 
co O 
P P 

o a 

(D « 
P (D 

S 


co 

o 


cc ^ & 

P co O 

i g s 


Proc. Linn. Soc. N.S.W., 132, 2011 



RELATIONSHIPS BETWEEN GEODIVERSITY AND VEGETATION 



Geological Substrates 



aeolian (red) sandplains 
estuarine sediments 
felsic intrusives 
felsic volcanics 
floodplain alluvium 
high quartz sedimentary 
lacustrine sediments 
laterite 



limestone 

low quartz sedimentary 

mafic volcanics & intrusives 

residual alluvial clay 

residual alluvial sands 

residual alluvial/colluvial sand & gravel 

siliceous (white) sandplains 

ultramafic igneous & metamorphics 


Figure 2. Geological map of New South Wales showing the distribution of 16 substrate types (see Table 
2 for description). 
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Climate surfaces 

Spatial grids for monthly precipitation, 
temperature and radiation parameters were generated 
by fitting spline functions (Hutchinson 1991) 
to weather station data across Australia. Using 
ANUCLIM v6.1 (http://fennerschool.anu.edu.au/ 
publications/so ftware/anuclim.php#overview), the 
weather station data were aggregated to monthly 
averages for a 30-year period (1975-2005) and 
interpolated across the continent on a 9-second 
latitudinal-longitudinal grid (approx. 250 m). The 
weekly climate averages were used to calculate the 
following parameters: MeanTemp- Mean temperature 
across all weeks of the year; MinTmp- mean of 
lowest weekly minimum temperature; MaxTmp- 
mean of highest weekly maximum temperature; 
DiRngTmp- mean of weekly diurnal temperature 
ranges; AnnRain- Mean annual rainfall; SummRain 
- Mean rainfall of December-February; RainDryMth- 
Mean rainfall of the driest month; MIdry- Mean 
moisture index of the driest month. Moisture index is 
calculated from weekly rainfall, evaporation and soil 
moisture storage (see http://fennerschool.anu.edu. 
au/publications/software/anuclim/doc/params.html). 
Spatial data for these parameters were re-projected 
to the Australian Geodetic Datum 66 and Lamberts 
Conformal projection in a geographic information 
system. 

Data analyses 

To examine the fidelity between vegetation types 
and geological substrates, 500 randomly located 
points were generated in a GIS to sample each of 
the 100 mapped vegetation classes. Subsequently, 
these were randomly sub-sampled to obtain 1000 
points in each of the 16 vegetation formations and 
subformations. The points sampling vegetation 
formations and classes were intersected with the 
geology map (i.e. vegetation cross-tabulated with 
geology) to estimate their frequencies of occurrence 
on each substrate type. Similarly, 1000 randomly 
located points were generated to sample each of the 
16 mapped geological substrates, and these were 
intersected with the vegetation map to estimate 
the diversity of vegetation formations and classes 
represented on each substrate type. Calculations 
of the number of points represented in the cross¬ 
tabulations were based on 90 lh percentile of points to 
reduce the effect of boundary errors in mapping that 
may cause spurious occurrences of vegetation types 
on particular substrates. Thus, the number of substrate 
types per vegetation unit was taken as the minimum 
number of substrates accounting for 90% of points 


within that unit. Pearson correlation coefficients were 
calculated to assess the association between the area 
of vegetation classes and formations and the number 
of substrates that they occupy. 

The variation in species composition among 
vegetation classes attributable to geological substrate 
and climate was evaluated using ordinations and 
partial variance analyses (Leps & Smilauer 2003). 
This required construction of three data matrices 
characterising the species composition, geological 
substrates and climatic habitat of vegetation classes. 
A presence/absence species matrix (99 classes x 
1625 species) was constructed from the floristic 
descriptions (vascular flora) of vegetation classes 
in Keith (2004), which had been compiled from 
frequently mentioned species and identified dominant 
species in the descriptions of source map units. A 
substrate matrix (99 classes x 16 geological substrates) 
was constructed from the relative representation of 
substrates types within each mapped vegetation class, 
estimated from frequencies of 500 random points 
per class, as described above. A climate matrix (99 
classes x climate variables) was constructed from 
mean climate parameters across the same samples 
of 500 random points, which were intersected with 
each of the nine climate surfaces. Unfortunately, 
climate surfaces were unavailable for offshore areas, 
so four vegetation classes (Oceanic Rainforests, 
Oceanic Cloud Forests, Coastal Headland Heaths and 
Maritime Grasslands) were omitted from the analysis. 
Hence the species matrix was reduced to 95 classes x 
1500 species. 

To determine whether species responses best 
fitted a linear or unimodal response to environmental 
gradients, redundancy analysis and canonical 
correspondence analyses were each carried out on the 
species matrix constrained by the combined matrices 
for geological substrates and climate parameters. The 
first four Eigen vectors of each analysis accounted 
for a similar proportion of floristic-environmental 
relationships (25%), although the redundancy analysis 
accounted for slightly more floristic variation than the 
canonical correspondence analysis (7.8% cf. 7.2%). 
A linear response model was therefore assumed for 
subsequent analyses. 

An unconstrained principal components 
analysis was first carried out on the species matrix. 
This allowed display of floristic relationships and 
examination of environmental relationships using 
indirect gradient analysis to fit vectors representing 
each variable in the substrate and climate matrices to 
the floristic ordination. 

Two partial redundancy analyses were then 
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carried out to quantify the proportion of Holistic 
variation uniquely attributable substrate types and 
climatic variables. These were done as constrained 
ordinations with substrate and climate defined as the 
environmental and covariablc matrices, respectively, 
and then vice versa for the second analysis (ter Braak 
& Smilaucr 1999). A third partial redundancy analysis 
was then carried out to determine the combined 
(union) proportion of floristic variation attributable 
to either substrate or climate with both sets of 
variables combined within a single environmental 
matrix. Finally, the (intersection) proportion of 
variation attributable to both substrate and climate 
combined was calculated by subtracting the sums of 
all canonical eigen values from the first two partial 
redundancy analyses from the sum of canonical eigen 
valyes in the third partial redundancy analysis (Leps 
& Smilauer 2003). All redundancy analyses were 
carried out in CANOCO for windows v4.02 (ter 
Braak & Smilauer 1999). 

RESULTS 

Fidelity of vegetation types to geological 
substrates 

While none of the vegetation formations was 
restricted to a single geological unit, there was a 
strong non-random relationship between vegetation 


formation and geological substrate (y 2 = 21562, 
P«0.001, df = 225). Approximately 27-61% of the 
extent of each formation occurred on one substrate. 
More than 90% of the distribution of ten of the 16 
formations was restricted to five or less geological 
substrates and all but one formation was restricted 
to seven or fewer substrates (Fig. 3). The number of 
substrates occupied by a vegetation formation was 
unrelated to the extent of its distribution (R—0.078, 
P>0.5, df=15), indicating that vegetation formations 
were not restricted to a small number of substrates 
simply by virtue of small distributions. Furthermore, 
geological substrates supported a limited range of 
vegetation formations. Four or fewer vegetation 
formations accounted for more than 90% of the area 
covered by eleven of the 16 substrate types and none 
of the substrate types supported more than seven 
formations. The number of formations per substrate 
type was unrelated to substrate area (R=0.43, P-0.1, 
df=15). 

At the level of vegetation class, the association 
with geological substrates was even more strongly 
expressed. Seven vegetation classes were essentially 
restricted to a single substrate type, while two-thirds 
of the 99 classes had at least 90% of their distribution 
restricted to three or fewer substrates (Fig. 4). As for 
formations, the number of substrates occupied by 
a vegetation class was unrelated to the extent of its 
distribution (R=0.027, P>0.5, df=98). 
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03 




Number of substrates per vegetation formation 


Figure 3. Fidelity of 16 vegetation formations to 16 geological substrate types. 
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Number of substrates per wgetation class 


Figure 4. Fidelity of 99 vegetation classes to 16 geological substrates. 


Rainforests, both Wet Sclerophyll Forest 
subformations and Grassy Woodlands were 
strongly associated with low-quartz sediments 
and metasediments (Appendix 1). Rainforests also 
occurred frequently on felsic volcanics, while the 
wet sclerophyll forests and grassy woodlands were 
more strongly associated with felsic intrusives. It 
is likely that the rainforests, grassy woodlands and 
grasslands were also well represented on mafic 
volcanics, but much of this substrate has been cleared 
of its native vegetation. The Heathlands and shrubby 
subformation of Dry Sclerophyll Forests were 
strongly associated with high-quartz sediments and 
siliceous (white) sands of marine origin, but also had 
significant representation on low-quartz sediments. 
In contrast, the shrub/grass subformation of Dry 
Sclerophyll Forests was primarily associated with 
low-quartz sediments and felsic intrusives (Appendix 
1). The Alpine Complex occurred mainly on felsic 
intrusives and low-quartz sediments. All three 
wetland formations and the Grasslands were strongly 
associated with active fluvial alluvium, with lower 
frequencies of occurrence across a range of other 
substrates. The shrubby Semi-arid Woodlands and 
Arid (acacia) Shrublands were strongly associated 
with aeolian (red) sands, while the grassy Semi-arid 
Woodlands occurred primarily on floodplain alluvium 
and residual alluvial clays and Arid (chenopod) 
Shrublands were on aeolian sands and residual clays 
(Appendix I). 

The vegetation classes that were essentially 
restricted (>90% of occurrence) to one substrate type 


included a range of Rainforests, Dry Sclerophyll 
Forests, Heathlands and Semi-arid Woodlands 
(Appendix 2). The geological substrates that supported 
the broadest ranges of vegetation formations include 
low- and high-quartz sedimentaries, felsic intrusives 
and floodplain alluvium (Appendix 1). 

Relative influence of geology and climate on 
vegetation 

The Principal Components ordination showed 
that vegetation classes within the same formation 
generally clustered together (Fig. 5a). This suggests 
considerable floristic affinities within formations, 
even though classes were grouped together within 
formations on the basis of structural and functional 
resemblance, rather than compositional resemblance. 
Indirect gradient analysis showed that geological 
substrates and climatic parameters account for a 
diverse array of compositional gradients within 
vegetation of the region (Fig. 5b). Individual climatic 
parameters appeared to exert a stronger influence 
on vegetation, as their vectors were generally 
longer than those representing individual geological 
substrates, indicating stronger correlations with 
species composition. However, geological substrates 
appeared to exert a more diverse range of influences, 
as their vectors spanned a greater range of directions 
than those representing climate parameters (Fig. 5b). 

Partial variance analysis showed that, in 
combination, the full set of geological substrates 
accounted for a greater proportion of variation in 
species composition than the climate variables (Fig. 6). 
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• RF 

• WSFs 
oWSFg 

- GW 
x GL 

a DSFg 
a DSFs 
+ HL 
x AC 

- FrW 

- FoW 

° sw 

O SAWg 

• SAWs 
■ ASc 

a ASa 



Figure 5. (a) Scatter plot of unconstrained Principal Components Analysis of 99 vegetation classes 
grouped by formations: RF- Rainforests, WSFs- Wet Sclerophyll Forests (shrubby subformation), WSFg- 
Wet Sclerophyll Forests (grassy subformation), GW- Grassy Woodlands, GL- Grasslands, DSFg- Dry 
Sclerophyll Forests (shrub/grass subformation), DSFs- Dry Sclerophyll Forests (shrubby subformation), 
HL- Heathlands, AC-Alpine Complex, FrW- Freshwater Wetlands, FoW- Forested Wetlands, SL- Saline 
Wetlands, SAWg- Semi-arid Woodlands (grassy subformation), SAWs- Semi-arid Woodlands (shrubby 
subformation), Arid Shrublands (chenopod subformation), Arid Shrublands (acacia subformation), (b) 
Plot of vectors representing 16 geological substrates (thin black lines) and 9 climate parameters (thick 
grey lines) fitted to the Principal Components ordination. The substrate types are: ResSand- Residual 
Alluvial Sand, Colluvial- Colluvial/alluvial sand and loam, ResClay- Residual alluvial clay, ActAlluv- 
Active alluvium, Lacust- Lacustrine sediments, EstSeds- Estuarine sediments, MarSand- Marine sands, 
Hqseds- High-quartz sedimentary rocks, Lqseds- Low-quartz sedimentary & metamorphic rocks, 
Felslnt- Felsic intrusives, FelsVolc- Felsic volcanics, UltMaf- Ultramafic volcanic and metamorphic 
rocks, Limestn- Limestone, MafVolc- Mafic volcanics. The climate parameters are: DiRngTmp- Diurnal 
range of temperature, MaxTmp- Mean temperature of the warmest month, MeanTemp- Mean annual 
temperature, MIdry- Mean moisture index (see text) of the driest quarter, MinTmp- Mean temperature 
of the coldest month, AnnRain- Mean annual rainfall, SummRain - Mean rainfall of December-Febru- 
ary, RainDryMth- Mean rainfall of the driest month. Note the vector plot (5b) is enlarged by a factor of 
2 relative to the scatter plot (5a). 
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Figure 6. Venn diagram showing portions of variation in floristic composition of vegeta¬ 
tion classes attributable to substrate only (S), climate alone (C), both substrate and climate 
(SC) and unexplained variation (U), as determined by partial redundancy analysis. 


Geological substrates and climate accounted for 
largely independent components of variation in 
species composition, as only 4.4% of total floristic 
variation was correlated with both geology and 
climate in combination. Together, geology, climate 
and their overlapping component accounted for just 
over one-third of total floristic variation, leaving two- 
thirds unexplained (Fig. 6). 

DISCUSSION 

The influence of geodiversity on vegetation 

Vegetation exhibited strong relationships with 
geodiversity at both class- and formation level 
across 80 million hectares of south-eastern Australia. 
Almost one-fifth of floristic variation across this 
large temperate region was uniquely attributable 
to geological substrates, independent of climatic 
variables. Each vegetation formation and class 
showed strong fidelity to a small range of geological 
substrates, with some classes restricted to a single 
substrate type. Stronger fidelity at the class level, 
relative to vegetation formations, indicates that 
relationships between vegetation and geodiversity 
are scale-dependent. At finer levels of vegetation 


classification than class, a still greater proportion of 
plant assemblages are restricted to a single type of 
substrate (e.g. Tozer et al. 2010). 

Indirect gradient analysis showed that species 
composition of vegetation was more strongly 
correlated with individual climate parameters, 
notably rainfall, than any single substrate type. 
However, the compositional trends associated with 
substrates encompassed a broader array of gradients 
than those associated with climate parameters. As a 
consequence, partial variance analysis showed that 
the substrate types collectively accounted for more 
variation than a set of parameters encompassing the 
means and extremes of climatic moisture, temperature 
and patterns in their seasonality. The overlapping 
component of floristic variation attributable to both 
geodiversity and climate was remarkably small. 

Amongtheclimatevariables,floristicrelationships 
with the three rainfall parameters were positively 
correlated with one another and negatively correlated 
with vectors representing maximum temperature, 
diurnal range and seasonality. This major gradient was 
associated with the transition from forested vegetation 
classes to semi-arid woodlands and arid shrublands. 
Minimum temperature and moisture index of the 
driest month (which incorporates evapo-transpiration 
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as well as precipitation) displayed somewhat different 
llorislic trends. 

A strong contrast was evident between substrates 
that produce impoverished soils (high-quartz 
sediments, marine (white) sands) and those that 
produce more fertile soils (low-quartz sediments, 
lelsic inlrusivcs and volcanics). The former were 
strongly associated with vegetation types dominated 
by sclcrophyllous shrubs (as understorey or canopy 
species), while the latter were associated with 
vegetation types with abundant mesophyllous 
shrubs and/or grasses. The mafic volcanic substrates 
generally define the upper limit of this soil fertility 
gradient, while vectors representing substrates with 
extreme levels of some mineral elements (ultramafics, 
limestone) are intermediate between those of mafic and 
felsic substrates. A similar but more subtle distinction 
is evident between vegetation classes found within 
dry-climate regions. Shrubby semi-arid woodlands 
and arid (acacia) shrublands are associated primarily 
with impoverished aeolian (red) sands and residual 
alluvial sands, clays and colluvium, while grassy 
semi-arid woodlands and arid (chenopod) shrublands 
are more common on active alluvium and residual 
clays. Estuarine and lacustrine sediments are uniquely 
associated with various types of wetlands, which are 
also associated with active fluvial alluvium. 

Support for biogeographic landscape theories 

The patterns described above are consistent with 
early comparative work between the flora of low- and 
high-quartz substrates in the Sydney region (Beadle 
1953, 1966) and with soil-vegetation relationships 
inferred from early survey work in western New 
South Wales (Beadle 1948). This work highlighted the 
association between sclerophylly and soil nutrients, 
notably phosphorus, which are more abundant in clay 
minerals derived from mafic substrates than felsic 
substrates and least abundant in quartz-rich substrates 
(Table 2). 

The observed vegetation-substrate patterns 
generally support Hopper’s (2009) characterisation of 
two general landscape types: Young Often Disturbed 
Fertile Landscapes (YODFELs) and Old Climatically 
Buffered Infertile Landscapes (OCBILs). ‘Young’ 
and ‘old’ in Hopper’s sense refer to age of landscape, 
rather than underlying geology. Hence YODFELs are 
characterised by relatively fertile soils whose nutrient 
capital has not been greatly depleted by leaching and 
which may undergo frequent disturbance related to 
fluvial or maritime events or mass movement. Their 
flora is dominated by recently evolved species with 
long-distance dispersal capabilities, propensity for 


colonisation, extensive distributions, generalist 
nutritional and reproductive biology, and tolerance 
of disturbance (Hopper 2009). The YODFEL 
profile fits many species of the grassy vegetation 
formations and subformations, which occur on the 
more fertile substrates (c.g. low-quartz sediments, 
volcanics, active alluvium). It also generally fits a 
large portion of the flora that characterises the three 
wetland formations, which may generally be viewed 
as occupying resource-rich sinks within regional 
landscapes (Keith 2004). 

In contrast, OCBILs are characterised by a 
diversity of ancestral species lineages with limited 
dispersal and colonisation capability, often with 
restricted distributions, specialised nutritional and 
reproductive biology, prominent sclerophylly and 
limited resilience to physical disturbance. Additional 
species traits associated with the sclerophyll syndrome 
were described in mechanistic detail and for a broader 
range of biota by Orions & Milweski (2007) in 
their “Nutrient-Poverty/Intense-Fire Theory”. The 
OCBIL profile describes many of the sclerophyll 
plant species that characterise substrates associated 
with impoverished soils (e.g. high-quartz sediments, 
leached marine sands, aeolian sands). Both landscape 
types appear to extend throughout the humid - arid 
climatic gradient of the region. 

It is noteworthy that much of rainforest flora 
does not readily fit either profile. Many of the taxa 
occupy climatically buffered environments and 
belong to ancient lineages that generally lack recent 
radiation and have suffered numerous extinctions 
(Crisp et al. 2004). Yet their habitats are not the most 
nutrient-impoverished nor very ancient landscapes 
and many of the taxa are widely dispersed with large 
distributions, some are ready colonisers. 

Axiomatic to both Nutrient-Poverty/Intense-Fire 
and OCBIL theories is the proposition that plants 
growing on nutrient-deficient soils with periodically 
adequate moisture, can synthesize ‘excessive’ 
carbohydrates, which are deployed to produce well- 
defended foliage, large quantities of lignified tissues 
and readily digestible exudates (Orions & Milewski 
2007, Hopper 2009). The nutritional properties of 
geological substrates therefore define a fundamental 
basis for evolution of Australian biota and retain a 
distinctive signature on the present-day distribution of 
vegetation formations and assemblages in the region 
of south-eastern Australia examined here. Given their 
strong influence on contemporary vegetation patterns, 
geological substrates which, with few exceptions, are 
essentially fixed landscape features over millenial 
time scales, appear to impose significant constraints 
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on vegetation response to climate change, especially 
in landscapes with OCBIL characteristics. 

Approximately two-thirds ofthefloristic variation 
remained unexplained in the direct gradient analysis. 
Part of this unexplained variation may include 
unrepresented influences of soils and climate. For 
example, substrate types were defined very broadly 
and often encompass considerable heterogeneity, not 
only in the complexes of rocks juxtaposed within 
them, but in the mineral composition and texture and 
structure of soils produced across catenary sequences 
of the landscape. The movement and availability of 
water across the landscape is also an important source 
of variation that is not fully represented by the climatic 
variables included in the current analysis. This 
essential resource almost certainly accounts for some 
of the unexplained floristic variation, particularly in 
the wetland component of the biota. 

Fire regimes are also likely to account for a 
fraction of the unexplained variation, as a lack of 
suitable spatial data precluded any consideration of 
them in the analysis. Fire regimes have been identified 
as driving evolutionary forces in Australia and other 
continents (Bond 2005, Bowman et al. 2009). They 
are an important component of Nutrient-Poverty/ 
Intense-Fire theory, as rapid accumulation of nutrient- 
poor biomass, a result of low rates of herbivory, 
provides fuel for intense fire, which in turn promotes 
nutrient poverty through volatilisation (Orions & 
Milewski 2007). Any remaining variation in floristic 
composition of south-eastern Australia is mostly 
attributable to sampling error and inherent spatial 
autocorrelation, as time lags in vegetation dynamics 
and limited dispersal processes impose an inherently 
clustered spatial structure on the composition of biota 
in the landscape. 

Map-based approach to ecological analysis 

The map-based approach employed in this 
study has both strengths and limitations. A major 
advantage is that it permits a balanced stratified 
random sampling design across the entire study area. 
This overcomes a significant constraint for analyses 
based on field samples over such a large region 
- the available data are inevitably skewed and non- 
randomly distributed across the landscape to varying 
degrees. A complementary analysis based on field 
samples may nonetheless be profitable, as it permits a 
more direct location-based exploration of vegetation- 
environment relationships, and hence exploration of 
finer-scale patterns than can be represented on sub¬ 
continental maps. 

A potential limitation of the map-based approach 
is that imprecision in the boundaries of both maps 


may have resulted in some combinations of vegetation 
and substrate types that do not occur in nature, as 
well as a margin of error in estimated frequencies 
of association. Non-concurrence of soil, soil parent 
material and bedrock could occur, for example, where 
there is significant lateral movement of sediment 
downslope from its origin. This promotes a tendency 
for the fidelity of vegetation types to substrate types 
to be under-estimated (i.e. vegetation types are more 
restricted to vegetation types than the data indicate). 
To offset such effects, frequencies in Figs. 3 and 4 
were based on the 90 lh percentile of sampled points, 
although it is uncertain whether this adjustment 
adequately compensated spatial errors in the absence 
of field validation data. A second limitation of the map- 
based approach is that it does not allow relationships 
between floristics and environmental data to be 
explored directly. This was because the species and 
substrate matrices were based on descriptive data 
averaged across the mapped range of each unit, rather 
than location-specific estimates. Thirdly, depending 
on the methods employed to generate source maps of 
for corresponding areas, the spatial data for vegetation 
and geology may not be independent throughout the 
mapped area. For example, in some cases geological 
boundaries may have been used as proxies for 
vegetation mapping and conversely remote sensing of 
vegetation may have been used to identify geological 
boundaries. Such non-independence may inflate map- 
based correlation between vegetation and geology. 
However, such effects are mitigated by the use of 
multi-criteria in remote sensing and modelling, only 
some of which will be non-independent proxies, as 
well as varying levels of field sampling to directly 
verify mapped units. Independent reclassification of 
the vegetation and geological maps further reduced 
any non-independence. While these methodological 
issues limited the resolution of relationships that could 
be examined, the analytical methods employed were 
sufficiently sensitive to detect major influences of 
geological substrates on vegetation that, collectively, 
appeared to be stronger than, and largely independent 
of climatic influences. 
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Appendix 1. Percentage of 1000 random points in each vegetation formation on each geological substrate type 
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Appendix 2. Percentage of 500 random points in each vegetation class on each geological sub¬ 
strate type (see Keith 2004 for description of classes). 
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